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Abstract 

The contribution of different components of an air shower to the total energy deposit 
in the atmosphere, for different angles and primary particles, was studied using the 
CORSIKA air shower simulation code. The amount of missing energy, parameter- 
ized in terms of the calorimetric energy, was calculated. The results show that this 
parameterization varies less than 1% with angle or observation level. The depen- 
dence with the primary mass is less than 5% and, with the high energy hadronic 
interaction model, less than 2%. The systematic error introduced by the use of just 
one parameterization of the missing energy correction function, for an equal mix- 
ture of proton and iron at 45°, was calculated to be below 3%. We estimate the 
statistical error due to shower-to-shower fluctuations to be about 1%. 
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1 Introduction 



In the last years, the quest to unravel the many mysteries related to the cosmic 
radiation has been intensified, in an attempt to answer open questions [1,2] 
about origin, propagation and chemical composition of the ultra high energy 
cosmic rays. The Auger [3] observatory, now under construction, and future 
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experiments like Telescope Array [4], EUSO [5] and OWL [6] will shed light on 
these questions. 



A common feature of these experiments is the use of the fluorescence technique, 
first explored by the Fly's Eye group [7] . The fluorescence telescopes use the 
atmosphere as a calorimeter, making a direct measurement of the longitudi- 
nal shower development, which represents the most appropriate technique to 
determine the energy of the primary particle. This is based on the assumption 
that the fluorescence yield is proportional to the local energy deposit. Recent 
measurements [8,9] have shown that this is valid for electrons and efforts are 
being made to improve [10,11,12] and extend the results. 

However, a fraction of the primary energy cannot be detected because it is 
carried away by neutrinos and high energy muons that hit the ground. Cor- 
rections for the so called missing energy must be properly applied to the 
calorimetric energy E ca \ measured, in order to find the primary energy Eq. 
The first parameterizations of the missing energy, as a function of the calori- 
metric energy, was done by J.Linsley [13,14] and the Fly's Eye group [15]. Some 
years ago, C. Song et al [16] recalculated it using simulated air showers. The 
method consisted of first determining the number of charged particles in the 
shower N c h(t), as a function of the atmospheric depth t, including a correction 
to take into account particles discarded below the simulation threshold. The 
track length integral was then evaluated, assuming a mean energy loss rate 
(a), to find an estimate of the total calorimetric energy: 



Nowadays, a very detailed simulation of the energy deposited in the atmo- 
sphere by air showers [17] is available in the Monte Carlo program COR- 
SIKA [18]. In this paper, we use it to study the contribution of different com- 
ponents of an air shower to the total energy deposit. This approach is different 
from the previous one [16] since we use directly the energy deposited in the 
atmosphere, by each component of the shower. 

With this information, we calculate the amount of missing energy and param- 
eterize the fraction E ca i/Eo as a function of E ca \. We derive this function for 
proton and iron primaries and compare it with previous results. Additionally, 
we investigate the dependence of this parameterization with the zenith angle, 
high energy hadronic interaction model and observation level. 

The outline of this paper is as follows: section 2 gives a brief description of the 
simulations performed. In section 3 we discuss the concepts of missing energy 
and energy deposit as ingredients of our calculations to obtain the primary 
energy. In section 4, we present a summary of the longitudinal profiles of 
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energy deposit generated in the simulations and we discuss the missing energy- 
correction curve obtained. Our final conclusions are presented in section 5. 



2 Simulation 

The CORSIKA simulation code, version 6.30, was used to study the energy 
balance in the propagation of atmospheric showers. This Monte Carlo package, 
designed to simulate extensive air showers up to the highest energies, is able 
to track in great detail almost all particles. 

In particular, this program provides tables for the energy dissipated in the at- 
mosphere as a function of the atmospheric depth [17]. For each kind of particle 
(gamma, electron, muon, hadron and neutrino), the energy lost by ionization 
in air and by simulation cuts (angle and energy thresholds) are specified. In 
addition, there is also information about the amount of energy carried by 
particles hitting the ground. 

In this paper, QGSJET01 [19,20] has been employed for higher energies E iab > 
SOGeV, and GHEISHA [21] for lower energies. Optimum [22,23] thinning [24] 
was applied, so as to reduce CPU time, while keeping the simulation detailed 
enough for our purposes. We have used a thinning level of Eth = 10" 4 , with 
maximum weight limitation set to u max = 10~ 4 E /GeV for muons, hadrons 
and electromagnetic particles. The energy threshold for electromagnetic par- 
ticles was set to E^ = 50keV and, for hadrons and muons, E\ hr = E^ hr = 
50MeV. Ground was fixed at sea level. 

Proton and iron initiated air showers were simulated for three different energies 
10 18 eV, 10 19 eV and !0 20 eV and four different angles 0°, 30°, 45° and 60°. For 
each combination of these parameters, 100 showers were simulated and the 
longitudinal profiles of energy deposit were carefully studied. 

With the angle fixed at 45°, we have also simulated proton and iron showers 
for two different conditions. First, with a different high energy hadronic inter- 
action model (SIBYLL2.1 [25,26,27]). Second, with a different detection level 
(300m above sea level). 



3 Energy Balance 

We will use a nomenclature for energy distribution that is similar to that usu- 
ally seen in analytic solutions of diffusion equations. The number of particles of 
type i at depth t, with energy in the interval [E, E + dE], is N^E, t)dE. There- 
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fore, Fi(E, t)dE = ENi(E, t)dE is the total energy of these particles. The inte- 
gral form of these distributions is written as F^E < E Q ,t)dt — Jq° F^E, t)dE 
and so on. 

For the shower processes, we define i = h (for hadrons), \i (for muons), v 
(for neutrinos), em (electromagnetic component). For the i t h component, the 
energy deposited due to ionization prior to depth t is given by Di(t). 

Then, at any depth, we can write an expression for the primary energy: 
E = F i (E>0,t)+ ]T Di(t) 

i=h,fi,v,em i=h,fi,v,em 

Now, if we introduce energy thresholds 1 we have 

E = £ F t (E > El hr , t) + £ F t (E < El hr , t) + £ D t (t) 

i i i 

For practical reasons, particles submitted to threshold cuts (both in energy 
and in angle) are no longer tracked, and this expression is not exact. CORSIKA 
computes the energy of these particles as dissipated energy. The exact equation 
must consider the sum of the energies of all particles below the threshold 
F*(E < El hr ,t), prior do depth t, and the energy deposited only by particles 
above the threshold D*{t): 

E = Y,F(E> El hr , t) + ^F*(E< El hr , + E A* W 

i i i 

Now we can define calorimetric and missing energy for our purposes. We take 
as calorimetric energy the energy deposited in air by ionization. The amount 
of energy which does not cause ionization is the missing energy. The energy in 
the electromagnetic component and a fraction f h of the energy in the hadronic 
component, no matter the stage of shower development, are also taken as 
calorimetric. This has to be so, because they are accounted for when evaluating 
the track length integral. 

Besides that, when a particle reaches a threshold, its energy is dropped from 
the simulation. However, a fraction of this energy would in reality go to ion- 
ization. We can take these fractions as constants, for each component, as 
they are averaged over a large number of low energy particles. Therefore, 
fiF*(E < El hr ,t) is the fraction of the sub-threshold particles energy which 

1 Even though there are also angular cuts, we write only E < E t h r , for simplicity. 
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goes to ionization. The remaining (1 — /$) goes to missing energy. The values 
of the constants and f' h will be given later. 

From the above equations, we can define the calorimetric energy at depth t 
as: 

Ecai(t) = F em (E > E thr , t) + F* m (E < E thr , t) + 
f' h F h (E > E? hr , t) + f h F*(E < E? hr , t) + 
f»F;(E<E? hr ,t) + J2D*(t) 



The remaining terms will add to the missing energy: 

E miss (t) = (1 - f h )F h (E > E h thr , t) + (1 - f h )F*(E < E? hr , t) + 
F,(E > E? hr ,t) + (1 - f,)F;(E < E? hr ,t) + 
F V {E > 0,t) 

We will consider now the fractions going to calorimetric and missing energies 
for particles submitted to threshold cuts. It is important to know the type of 
particles contributing and since this information was not available, we modified 
the CORSIKA code. For each particle dropped from the simulation, its type, 
kinetic energy and releasable energy were recorded. The releasable energy is 
the energy that would have been lost had the particle not been dropped from 
the simulation. 

Thus we could calculate the contributions of the different types of particles 
and their respective energy spectrum. Table 1 shows the mean kinetic energy 
and the relative contributions calculated from 500 proton showers with 10 19 eV 
and 45°. The relative contribution is the ratio of the energy deposited by a 
given type of particle to the total energy deposited by the respective shower 
component. The thinning level was set to 10~ 4 without weight limitation. We 
found similar results for iron induced air shower simulations. 

This information was used to further track these sub-threshold particles using 
the Geant4 code [28]. For each type of particle, 5000 air shower simulations 
initiated by these low energy particles were performed. The energy spectrum 
obtained from CORSIKA was used to set the energy of these events and the 
secondary particles were followed through a sea level standard atmosphere 
until IkeV. The Geant package allowed us to find out the amount of energy 
lost by ionization and other processes, and determine the fractions of energy 
deposited by ionization. These fractions are also shown in table 1. For the 
low energy electromagnetic particles, we made the assumption that all the 
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energy is deposited (f em = 1). For muons and hadrons, we used the fractions 
fn = 0.425 and f h = 0.739. 

It is also necessary to consider what happens to the hadrons hitting the ground. 
They have a significantly different composition and energy spectrum from 
those hadrons cut in air and we cannot assume that f' h = fh- We used the 
CORSIKA information about particles at detection level to analyze the energy 
spectrum and relative contributions of different types of hadrons. Table 2 
shows the results for the same set of 500 CORSIKA showers. In this case, the 
average energies are higher and we found the dependence of f' h on energy to 
be a second order effect. Geant4 was then used to simulate 500 showers, for 
each particle, with energy fixed to its average kinetic energy. We found that 
fk = 0-61. 

Therefore, the prescription used for calculating the missing energy at ground 
level is: 

E miss = 0.390F h (E > E? hr ,grd) + 0.2QIF*(E < E* hr ,grd) + 

1.000F M (£ > E? hr , grd) + 0.575F;(£ < E? hr , grd) + (1) 

1.000F„(£ > 0,grd) 

Since E ca i/E = 1 — E miss /E , we can calculate the fraction of the primary 
energy measured with the fluorescence telescopes. 



4 Results 

For the sake of clarity, we present in this section the results from our sim- 
ulations. In table 3, we show the contribution for the energy deposit from 
different shower components, as a percentage of the primary energy, for each 
simulated angle and energy. 

Each three lines, corresponding to a fixed choice of angle and energy, refer to 
three different contributions to the energy deposit, which we called: (ION) ion- 
ization in air, (CUT) simulation cuts, and (GRD) particles arriving at ground. 
Each column gives the contribution from different shower components: gam- 
mas, electrons, muons, hadrons and neutrinos. The first value (first/second) 
refers to proton primaries, while the second refers to iron primaries. 

The first line accounts for ionization energy losses in air, consequently we have 
no contributions from gammas or neutrinos. We considered all contributions 
from this line as calorimetric energy. 
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The second line gives the total energy of particles left out of the simulations be- 
cause of angular or energy cuts. A fraction of this energy will go to ionization, 
as described in last section. Since neutrinos are not tracked in the simulation, 
they are dropped as soon as produced and all their energy appears in this line. 
We considered all their energy as missing energy. 

The last line gives the amount of energy carried by particles arriving at the ob- 
servation level. In this case, we considered all energy from the electromagnetic 
and a fraction of the hadronic component hitting the ground as calorimetric 
energy. The energy carried by muons was taken as missing energy. 



4-1 Missing energy 

We have used the data in table 3 to evaluate the mean missing energy for each 
combination of angle, energy and primary particle as given by equation 1. 

In figure 1, the correction factor for the missing energy is shown, plotted as 
the fraction E ca i/E Q , as a function of E cai . We can see that the variation in the 
amount of missing energy is only slightly dependent on the angle and decreases 
with energy. At lO 18 ^, the ratio E ca i/E varies between 90.8 — 91.3% for 
protons, and between 86.4 — 87.0% for iron. 

On the other hand, the ratio E ca i/E is largely dependent on the primary mass. 
At a fixed angle, the difference in the amount of missing energy for proton 
and iron showers decreases with energy. For instance, at 45° this difference 2 
is 4.9% at 10 18 eV, 3.3% at 10 19 eV and 2.1% at 10 20 eV. 

A mean parameterization of E ca i/E , taken as a mixture of 50% proton and 
50% iron at 45°, is usually used in the reconstruction routines. In table 4, we 
give the fitting parameters of the mean missing energy correction function for 
different simulations conditions. For showers with zenith angle 45°, we have: 



The systematic error in the value of E ca i/E , when calculating the primary 
energy with the mean parameterization, is less than ±2.5% ± 0.4% at lEeV 
and decreases with energy. Shower-to-shower fluctuations were evaluated, and 
found to be independent of energy. The rms value of E ca i is 1.1% of the primary 
energy for proton showers, and 0.4% for iron showers. 



2 In this and in the next section, all differences are relative, e.g., (v\ — V2)/v. 
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Figure 2 shows in detail the mean missing energy correction, calculated for a 
mixture of 50% proton and 50% iron. As can be noticed, there is a small shift in 
the curves relative to the previous results by Song [16] and Linsley [14]. Song's 
results are based on Monte Carlo simulations, and the mixed composition is 
shown. Linsley's results were calculated based on muon size measurements, 
and are independent of the primary composition. Results agree to within 1%, 
but they were not obtained under the same conditions. In figure 3, we compare, 
for 45°, the missing energy correction curves obtained for the ground level set 
to 0m and 300m, which is the same used by Song et al. As can be seen, our 
result does not depend on the observation level. 

The dependence on the interaction model is considerable, as shown in figure 
4. We compared, for proton and iron at 45°, the missing energy correction 
curves obtained with QGSJET01 and SIBYLL2.1 models. The difference is 
practically constant: around 1.6% for proton and 1.2% for iron. In all cases, 
SIBYLL predicts less missing energy than QGSJET. 



5 Conclusions 



The energy deposit by atmospheric air showers was studied aiming for a bet- 
ter reconstruction of the primary energy. The new and very detailed energy 
balance now present in CORSIKA was used. 

Our results for the missing energy correction function agree with previous 
calculations [14,16], to within 1%. The dependence of this function on angle 
and observation level was found to be less than 0.7%. Comparing the high 
energy hadronic interaction models QGSJET and SIBYLL, the results differ 
by less than 1.6%. The larger dependence comes from the primary mass, being 
less than 5% at 10 18 eV, decreasing with energy. 

We found a mean parameterization of E ca i/E , taken as a mixture of 50% 
proton and 50% iron at 45°, as usually used in the reconstruction routines. 
Considering the QGSJET model only, we estimate that the total systematic 
error, introduced by the use of just this parameterization, is below 3% at 
lEeV, and below 2% at lOOEeV. This is of the order of shower-to-shower 
fluctuations. For proton showers, the rms value of E ca i is 1.1% of the primary 
energy. For iron showers, it is 0.4%. 
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41.4 
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Table 1 

Particles discarded in air. Mean kinetic energy, relative contribution and fraction of 
energy going to ionization are shown. The fractions f em , and fh are the average 
over all particles weighted by their relative contributions. 



particle E kin (GeV) Rel. Contri. (%) Ion. Frac(%) 



n 


10. 


7.12 


70.1 ± 0.1 


V 


31.6 


5.57 


75.3 ± 0.1 


P 


100. 


3.21 


73.2 ± 0.1 


n 


100. 


2.92 


72.2 ± 0.1 


Kl 


1000. 


7.30 


36.3 ± 0.2 


K± 


1000. 


8.86 


60.4 ± 0.2 




316. 


32.06 


61.7 ± 0.2 


7T+ 


316. 


33.01 


59.3 ± 0.2 



fh = 61-0 



Table 2 

Hadrons at ground level. Mean kinetic energy, relative contribution and fraction of 
energy going to ionization are shown. The fraction f' h is the average over all particles 
weighted by their relative contributions. 
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gammas electrons muons hadrons neutrinos 







ION 


-/- 


64.7 / 65.3 


1.0 / 1.5 


0.2 / 0.3 


-/- 




lEeV 


CUT 


1.1 / 1.2 


10.7 / 10.8 


0.1 / 0.1 


0.3 / 0.5 


3.0 / 4.6 






GRD 


8.4 / 4.6 


4.2 / 2.1 


5.2 / 8.0 


1.1 / 1.1 


-/- 






ION 


-/- 


62.0 / 64.9 


0.8 / 1.2 


0.2 / 0.3 


-/- 


0° 


lOEeV 


CUT 


1.1 / 1.2 


10.2 / 10.7 


0.1 / 0.1 


0.3 / 0.4 


2.5 / 3.6 






GRD 


11.3 / 7.0 


5.9 / 3.4 


4.4 / 6.2 


1.1 / 1.2 


-/- 






ION 


-/- 


57.7 / 62.8 


0.7 / 1.0 


0.2 / 0.2 


-/- 




lOOEeV 


CUT 


1.0 / 1.1 


9.4 / 10.3 


0.1 / 0.1 


0.3 / 0.3 


2.1 / 2.8 






GRD 


15.2 / 10.1 


8.4 / 5.1 


3.7 / 4.9 


1.2 / 1.3 


-/- 






ION 


-/- 


69.8 / 67.6 


1.2 / 1.8 


0.2 / 0.3 


-/- 




lEeV 


CUT 


1.3 / 1.3 


13.0 / 12.7 


0.1 / 0.1 


0.4 / 0.5 


3.3 / 5.0 






GRD 


3.4 / 1.7 


1.5 / 0.7 


5.3 / 8.0 


0.5 / 0.4 


-/- 






ION 


-/- 


69.0 / 69.1 


1.0 / 1.4 


0.2 / 0.3 


-/- 


30° 


lOEeV 


CUT 


1.3 / 1.3 


12.9 / 12.9 


0.1 / 0.1 


0.4 / 0.4 


2.7 / 3.9 






GRD 


5.0 / 2.7 


2.4 / 1.2 


4.4 / 6.2 


0.5 / 0.5 


-/- 






ION 


-/- 


67.6 / 69.5 


0.8 / 1.1 


0.2 / 0.2 


-/- 




lOOEeV 


CUT 


1.3 / 1.3 


12.6 / 13.0 


0.1 / 0.1 


0.3 / 0.4 


2.3 / 3.1 






GRD 


7.0 / 4.0 


3.6 / 1.8 


3.7 / 5.0 


0.5 / 0.5 


-/- 






ION 


-/- 


71.2 / 67.2 


1.4 / 2.0 


0.2 / 0.3 


-/- 




lEeV 


CUT 


1.5 / 1.4 


15.6 / 14.7 


0.1 / 0.1 


0.4 / 0.5 


3.6 / 5.4 






GRD 


0.5 / 0.3 


0.2 / 0.1 


5.2 / 7.9 


0.1 / 0.1 


-/- 






ION 


-/- 


72.3 / 69.8 


1.2 / 1.6 


0.2 / 0.2 


-/- 


45° 


lOEeV 


CUT 


1.5 / 1.5 


15.8 / 15.3 


0.1 / 0.1 


0.4 / 0.5 


3.0 / 4.2 






GRD 


0.8 / 0.4 


0.3 / 0.2 


4.3 / 6.1 


0.1 / 0.1 


-/- 






ION 


-/- 


72.7 / 71.5 


1.0 / 1.3 


0.2 / 0.2 


-/- 




lOOEeV 


CUT 


1.5 / 1.5 


15.9 / 15.7 


0.1 / 0.1 


0.4 / 0.4 


2.6 / 3.4 






GRD 


1.2 / 0.7 


0.5 / 0.3 


3.7 / 4.8 


0.1 / 0.1 


-/- 






ION 


-/- 


67.5 / 63.4 


1.7 / 2.3 


0.2 / 0.2 


-/- 




lEeV 


CUT 


1.7 / 1.6 


19.2 / 18.1 


0.1 / 0.1 


0.5 / 0.6 


4.1 / 6.1 






GRD 


0.0 / 0.0 


0.0 / 0.0 


5.1 / 7.6 


0.0 / 0.0 


-/- 






ION 


-/- 


69.3 / 66.2 


1.4 / 1.9 


0.1 / 0.2 


-/- 


60° 


lOEeV 


CUT 


1.7 / 1.7 


19.7 / 18.8 


0.1 / 0.1 


0.4 / 0.5 


3.3 / 4.7 






GRD 


0.0 / 0.0 


0.0 / 0.0 


4.0 / 5.8 


0.0 / 0.0 


-/- 






ION 


-/- 


70.1 / 68.2 


1.2 / 1.6 


0.1 / 0.2 


-/- 




lOOEeV 


CUT 


1.8 / 1.7 


19.9 / 19.4 


0.1 / 0.1 


0.4 / 0.5 


2.9 / 3.8 






GRD 


0.0 / 0.0 


0.0 / 0.0 


3.5 / 4.6 


0.0 / 0.0 


-/- 



Table 3 



Mean energy deposit contributions (in percentage of primary's energy) from different 
shower components. Each three lines correspond to a fixed choice of angle and 
energy, and discriminate: (ION) ionization in air, (CUT) simulation cuts and (GRD) 
particles arriving at ground. The first value (first/second) refers to proton showers 
while the second refers to iron showers. 
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Fig. 1. Missing energy correction plotted as the fraction E cci i/Eq as a function of 
E ca i. The curves are for proton and iron showers at four different angles. 



Coefficients of correction function 







Iron 




Iron/Proton 




Proton 




angle 


A 


B 


C 


A 


B 


C 


A 


B 


C 





0.970 


0.100 


0.139 


0.977 


0.085 


0.117 


0.984 


0.071 


0.089 


30 


0.971 


0.102 


0.137 


0.979 


0.088 


0.114 


0.986 


0.074 


0.088 


45 


0.977 


0.109 


0.130 


0.967 


0.078 


0.140 


0.958 


0.048 


0.162 


60 


0.962 


0.098 


0.161 


0.948 


0.062 


0.220 


0.942 


0.035 


0.337 



Table 4 



Fitting parameters for different simulation conditions, as plotted in figure 1. The 
mid column indicates the 50%/50% mixture and the values corresponding to 45° 
are bold faced. The fit function used was: A - B(E/EeV)~ c . 
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Fig. 2. Missing energy correction plotted as the fraction E cci i/Eq as a function of 
E ca i. The average between iron and proton missing energy correction is plotted. For 
comparison, Song's and Linsley's results are also shown. 




Fig. 3. Missing energy correction plotted as the fraction E ca i/Eo as a function of 
E^. The dependence of the parameterization on ground level on ground level is 
shown. Simulations for proton and iron primaries, at 45°. 
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Fig. 4. Missing energy correction plotted as the fraction E ca i/Eo as a function 
of E ca i. The variation with the high energy hadronic interaction model is shown. 
Simulations for proton and iron primaries, at 45°. 
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